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In classical estimation theory, the central limit theorem implies that the statistical error in a 
measurement outcome can be reduced by an amount proportional to n _1//2 by repeating the measures 
n times and then averaging. Using quantum effects, such as entanglement, it is often possible to do 
better, decreasing the error by an amount proportional to n" 1 . Quantum metrology is the study 
of those quantum techniques that allow one to gain advantages over purely classical approaches. 
In this review, we analyze some of the most promising recent developments in this research field. 
Specifically, we deal with the developments of the theory and point out some of the new experiments. 
Then we look at one of the main new trends of the field, the analysis of how the theory must take 
into account the presence of noise and experimental imperfections. 



Any measurement consists in three parts: the prepa- 
ration of a probe, its interaction with the system to be 
measured, and the probe readout. This process is often 
plagued by statistical or systematic errors. The source 
of the former can be accidental (e.g. deriving from an 
insufficient control of the probes or of the measured sys- 
tem) or fundamental (e.g. deriving from the Heisenberg 
uncertainty relations). Whatever their origin, we can re- 
duce their effect by repeating the measurement and av- 
eraging the resulting outcomes. This is a consequence of 
the central-limit theorem: given a large number n of in- 
dependent measurement results (each having a standard 
deviation Act), their average will converge to a Gaussian 
distribution with standard deviation Aa/^/n, so that the 
error scales as n -1 / 2 . In quantum mechanics this be- 
havior is referred to as "standard quantum limit" (SQL) 
and is associated with procedures which do not fully ex- 
ploit the quantum nature of the system under investi- 
gation 1 . Notably it is possible to do better when one 
employs quantum effects, such as entanglement among 
the probing devices employed for the measurements, e.g. 
see Refs. [lHff. Consequently, the SQL is not a funda- 
mental quantum mechanical bound as it can be surpassed 
by using "non-classical" strategies. Nonetheless, through 
Heisenberg-like uncertainty relations, quantum mechan- 
ics still sets ultimate limits in precision which are typi- 
cally referred to as "Heisenberg bounds" . In Fig. [T] we 
present a simple example which may be useful to un- 



1 In quantum optics, the n~ 1 ' 2 scaling is also indicated as "shot 
noise" , since it is connected to the discrete nature of the radiation 
that can be heard as "shots" in a photon counter operating in 
Geiger mode. 



derstand the quantum enhancement. Part of the emerg- 
ing field of quantum technology [7], quantum metrology 
studies these bounds and the (quantum) strategies which 
allows us to attain them. More generally it deals with 
measurement and discrimination procedures that receive 
some kind of enhancement (in precision, efficiency, sim- 
plicity of implementation, etc.) through the use of quan- 
tum effects. 

This paper aims to review some of the most recent de- 
velopment of the field. For a more historical perspective 
of quantum metrology, we refer the reader to [5|. We 
start by introducing some important results on quantum 
estimation theory which focus on the optimization of the 
probe readout. Then we report some recent findings ob- 
tained in the context of parameter estimation for chan- 
nels, which allows also for the optimization of the probe 
preparation. It is only at this stage that Heisenberg-like 
scaling is obtained. Schemes based on filtering protocols 
and nonlinear effects will be also introduced. In the re- 
maining of the paper we will then deal with the analysis 
of the typically very fragile quantum metrology protocols 
in the presence of noise. 



Basics on Quantum Estimation for states 

In its simplest version a typical quantum estimation 
problem [a-[l4| consists in recovering the value of a con- 
tinuous parameter x (say the phase if of Fig. [1} which 
is encoded into a fixed set of states p x of a quantum 
system S. As in the example of Fig. [1] we can de- 
scribe it as a two step process where we first perform 
a measurement on S, and then extrapolate the value 
of x with some data-processing of the measurement re- 
sults. The measurement is described by a Positive Op- 
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FIG. 1: Ramsey interferometry. Its aim is to measure an unknown relative phase ip which is picked up by two orthogonal states 
| a), 1 6} of an atomic probing system. This is a generic framework that encompasses many different interferometric measurements 
(e.g. frequency-standards, magnetometry, optical phase, etc.) [H]. Left: In a conventional setup the probe preparation consists 
in producing each atom in the superposition |V>m) = (|a) + |6))/v2, which yields the output state \if) v ) = (|o) + e l¥> |6))/\/2 after 
the probing stage (black-boxes in the figure). The readout consists in checking whether \4> v ) is still in the initial state \ipi n }. 
This happens with probability p — \{ipi n \ip^}\ 2 = (1 — cosip)/2. Therefore by taking the ratio between the number of successes 
and the total number of readouts, we can recover the phase ip. If we repeat this measurement n times, the associated error on 
our estimation of ip, can then be evaluated by the standard deviation on the determination of p (i.e. Ap = \/p(l — p)/n) and by 

using error propagation theory, obtaining the SQL scaling Stp n = \J /\&p\ ~ n ~ 1 ^ 2 ■ Right: Quantum enhanced setting: A 
simple quantum strategy consists in dividing the n probes into groups of N, prepared in an entangled state (\a)® N + 16)® )/ y/2. 
Since each of the N vectors |6) acquires a relative phase ip, the final state is (\a)® N + e v \b) )/y/2. The probability that 



this state is equal to the initial one is now p e „t = (1 
can repeat this procedure v times, obtaining an error Sip n 



cos Nip) /2. Furthermore since we have v = n/N groups of probes, we 
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precision with respect to previous case, namely, the Heisenberg bound for phase sensing 0, E| 



(nN) 1//2 , with a N 1 ^ 2 enhancement in 



erator Valued Measurement (POVM) of elements 

{E^ n) } (where n is the number of copies of p x we use). 
The conditional probability of getting the outcome y 
will be then computed as p n (y\x) = Tr[E^ pf n ]. The 
data-processing of the result y will yield our estimate 
z of the value of x and, in the most general case, will 
be characterized by assigning some conditional proba- 
bilities p£st( z \y)- Ideally, we would like to have z as 
close as possible to the parameter x; unfortunately, de- 
pending on the physics of the problem and, possibly, 
on the selected estimation strategy, there will often be 
some residual uncertainty in the determination of the pa- 
rameter. This uncertainty is fully characterized by the 

probability P„(z\x) := J2 V Pest{ z \v)Pn{y\%), which de- 
scribes the statistical dependence of z on the true value 
x. It allows us to define the Root Mean Square Error 
(RMSE) SX n := P n (z\x) which pro- 

vides a good measure of the estimation accuracy (differ- 
ent measures are possible but will not be treated here, 
see e.g. Ref. [HI]). Under reasonable assumptions on the 
asymptotic behavior of the estimation function (i.e. re- 
quiring it to be asymptotically Locally Unbiased?), SX n 



2 This means that for n — > oo there exists a value of the param- 
eter x for which the estimator tends to the correct value, and 
at the same time the derivative of the estimator at that point is 
unity [To| . All reasonable estimators satisfy this condition while 
the pathological ones, do not (e.g. trivial estimation procedures 



can be shown to obey to the so called Cramer-Rao (CR) 
bound [HI], which implies 



SX n > l/y/F n {x), 



(1) 



where F„(x) := J2 y i^M^) /Pn(y\x) is the Fisher in- 
formation associated with the selected POVM measure- 
ment. 

Optimizing Eq. ((T|) with respect to all possible POVMs 
one then gets the inequality 



6X n ^ 



1 



1 



^/max £( „) F n (x) y/n J(p x ) 



(2) 



The term on the right is the quantum Cramer- Rao (q-CR) 
bound. It is obtained by exploiting an upper bound of 



that obtains the correct value of x by chance, as the case of a 
stopped clock which is correct twice per day). A generalization 
of Eq. {T} which applies to all estimation strategies without the 
assumption of locally unbiasedness was provided in 8, 9]. It is 
obtained by replacing 5X n on the lhs of Eq. |(TJ with its reg- 
ularized version SX n = \/y^. r \z(x) — x] 2 P n (z\x) in which z(x) 
is the quantity z divided by -^"^2 z z P n (z\x). Even though 
the inequality for 5X n refers to a larger number of estimation 
strategies, it is weaker since SX n ^ SX n . Note also that the 
regularization implies that SX n will diverge for the pathological 
estimators that would be excluded by the restrictions adopted in 
deriving JTJ|. 



3 



max £ („) F n (x) in terms of the quantum Fisher informa- 
tion 3 J(px). For instance, consider the case in which p x 
are pure states of the form \ip x ) = exp[— iHx] |0) with H 
an Hermitian operator and |0) a reference vector. Then, 
the q-CR bound takes the simple form of an uncertainty 
relation 



5X n ^ 



(3) 



where AH := y/((H - (H)) 2 } is the spread of H on |0). 

Three things are worth stressing at this point: i) the 
bound in Eq. @ holds for all possible POVMs, including 
those which operate jointly on the n copies while exploit- 
ing entanglement resources; ii) the SQL scaling n~ x / 2 
on the rhs is a direct consequence of the additivity of 
the quantum Fisher information when applied to tensor 
states pf n , i.e. J(pf n ) = nJ(p x ); Hi) in the asymptotic 
limit of large n, the q-CR bound is always achievable 
and the estimation strategy which attains it can be con- 
structed via local measurements and adaptive estimators, 
a strategy that employs only Local Operations and Clas- 
sical Communication (LOCC) Gj, [Mill ■ This im- 
plies that entangled resources at the measurement stage 
are not necessary to achieve the q-CR bound. Local mea- 
sures and some clever classical data-processing are suffi- 
cient. It also shows that the quantity 1/ y/nJ(x) has a 
clear operational meaning and can be used to quantify 
how hard the estimation problem is. 

Finally, we briefly consider the multi-parameter case 
where x is a vector of random variables |U>ll4l |. Also in 
this case one can construct an inequality in terms of the 
associated quantum Fisher information, but in general 
it will not coincide with the ultimate achievable bound. 
This inequality applies only to all n-body separable mea- 
surements [20l | - i.e. POVMs which may act globally on 
pf n but whose elements can be expressed as convex com- 
bination of tensor products of positive operators (which 
is a larger class than LOCC [21]). If, instead, we allow 
for generic joint measurements which exploit entangled 
resources in a non trivial fashion, then the bound may not 
hold and better performances are expected (even though 
the n _1//2 scaling still holds). 



Quantum parameter estimation for channels 

The analysis of the previous section focused on scenar- 
ios in which the set of states pf n is fixed by the estimation 
problem. Thus, it can only describe the last stage (probe 
readout) of a general measurement. Since (0) and ([3]) 



obey a SQL scaling n^ 1 / 2 , to obtain a quantum metrol- 
ogy type enhancement, we need to consider the other 
two stages (i.e. probe preparation and interaction). For 
this purpose it is useful to describe the correspondence 
x — > p x in terms of a quantum channel $ x which pro- 
duces p x when acting on the initial input state po of the 
probe, through the mapping $ x (po) = Px- Aim of this 
section is to characterize the best estimate of x which 
can be obtained for a fixed number n of applications of 
$ x , while optimizing with respect to the measurement, 
the estimation functions, and the choice of the initial 
state [22l - [3"7l | . This captures the basic aspects of most 
quantum metrology applications [f|. 
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FIG. 2: Schematic representation of parallel estimation 
strategies. In the picture, a white wedge element represents 
an input probe entering the apparatus, a black box represents 
an application of the channel <& x , a blue semicircles describes a 
local measurement on the probe, while a green box represents 
an entangling operation among the probes. CC (classical- 
input, classical-output): these are those strategies where the 
probes are prepared in a separable state (not necessarily 
all probes must be in the same state as depicted), and LOCC 
measurements at the output. In this case, the averaging of 
the local results yields a decrease in the result's precision that 
at most scales as the SQL, i.e. n _1//2 . CQ (classical-input, 
quantum-output): entanglement among the probes is gener- 
ated just before the detection. QC (quantum-input, classical- 
output): entanglement among the probes is generated before 
they are fed into the channel while no entanglement resource 
is employed at the detection stage which uses a LOCC strat- 
egy- QQ (quantum-input, quantum-output): the most gen- 
eral strategy where entanglement can be used both at the 
probe preparation and at the probe detection stages. By con- 
struction the QQ and CC strategies will always provide the 
best and the worst performances respectively. In the case of 
estimation of unitary channels, CQ has the same n -1 ^ 2 yield 
of CC, whereas QC and QQ can achieve the Heisenberg bound 
n _1 [2J]. Non-unitary channels have equal or worse perfor- 
mances (since any non-unitary map can be purified into a 
unitary), depending on their action on the probes [24l. I25I]. 



3 This is defined as J(p x ) := Tr[K^{p x )p x 1lp}(p' x )}, where p' x = 
dp x /dx, and where Kp 1 ^) := J2j,k:X j +\ k ^0 2O jk\j)( k \/{ x j + 
Afc) is the symmetric logarithmic derivative written in the basis 
that diagonalizes p x = £\ Xj\j)(J\ Ho| . 



As shown in Figs. [2] and [3l we can employ different 
strategies depending on the choice for the initial state 
of the probes, for the measurements, and for the oper- 
ation we can perform on the probed system (in parallel 
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or sequentially) |24[. In full analogy to what is done in 
quantum channel communication theory, we can also con- 
sider entanglement assisted schemes in which the incom- 
ing probes are entangled with an external ancillary sys- 
tem A which is not effected by the channel $ x (this is the 
most general method of performing measurements). Ac- 
curacy enhancements via entanglement assisted schemes 
for parallel configurations have been derived for special 
classes of quantum channels $> x [281432T ] . 
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FIG. 3: Right: Sequential (or multi-round) scheme [2^ - 
l26l . l38t |39( |. where a single probe state samples the n black 
boxes sequentially. Left: Entanglement-assisted sequential 
scheme, where an external ancillary system A aids the esti- 
mation. These schemes have some advantages with respect to 
the parallel schemes as in principle do not require to use multi- 
partite entanglement among n different probes to achieve sub- 
shot-noise scaling (furthermore they can simulate any other 
schemes by properly tuning the intermediate controls). The 
most general estimation scheme will be a combination of the 
sequential and of the parallel schemes of Fig. [2] An optimiza- 
tion over the possible black-box dispositions was presented 
in [33|. 

Optimizing with respect to all available resources can 
be a hard problem. For the sake of simplicity in the fol- 
lowing we will focus on the parallel strategies. In this 
context the conventional approach is to split the opti- 
mization into two stages [281 - 4321 ] . Specifically one first 

fixes a (possibly entangled) input state p^ and opti- 
mizes with respect to the measurements, and then mini- 

(n) 

mizes the result with respect to the p . From the results 
of the previous section it follows that the accuracy of the 
first stage will not be larger than q-CR bound associated 
with the selected state. Therefore, the second optimiza- 
tion yields 



SX n ^ min 



I 



J{P ( * ] ) 



(4) 



where J(p x 1 ^) is the quantum Fisher information of the 

output state p x n ^ . Since the quantity </(•) on the right is 
convex [28j . the minimum is always achieved by choos- 
ing pure states. In addition, due to the additivity prop- 
erties of </(•) it is clear that restricting the minimiza- 
tion in Eq. (j4]) to input states which are separable with 
respect to n probes, one can achieve at most a n -1 / 2 
scaling (SQL) (24|. Better scalings require entangled 

(n) 

states px . To exemplify this, we focus on the special 
case in which the & x perform random unitary rotations 



e~ zxH . This is a paradigmatic model which, among oth- 
ers, includes the Ramsey setup of Fig. [TJ and the Mach- 
Zehnder interferometer [3, [j| 0, E| , commonly used to 
detect an unknown phase shift (see below). Assume that 
the total number n — vN of probes is split in v 3> 1 
groups of N elements, where each group is prepared in 
the same entangled input state \ipQ N ^}. The output state 
of all n probes after passing through the black-boxes is 

:= (e-^*)®"!^)®". The bound i 
then be evaluated exploiting Eq. ^ and yields [24 



can 



5X n > 



(5) 



where I/c is the largest gap in the spectrum of the 
generator H of the transformation $ x . For N > 1 
it presents a A^-enhancement with respect to the 
standard-quantum-limit. The term on the right is the 
Heisenberg bound and can be attained by preparing the 

probes in each group into the entangled state \ipQ } = 

{\X max )® N +\X m in)® N )/V2, where |A max ), |A mm ) are the 
eigenvectors associated to the maximum and minimum 
eigenvalues of H, respectively 4 . From the previous sec- 
tion it follows that the threshold is asymptotically 
attainable for v 1 by exploiting POVMs which act lo- 
cally on the blocks and a maximum likelihood approach. 
Alternatively 24{ it is also achievable via adaptive strate- 
gies based on POVM which acts locally on each probe (so 
that the yield of QC strategies coincides with the one of 
QQ strategies). 

In contrast, the situation is more complicated for fi- 
nite v, specifically for v = 1 where Eq. ^ would 
yield the impressive n _1 scaling for the RMSE of the 
Heisenberg bound. For the special case in which $ x 
induces unitary transformations, this problem was an- 
alyzed in Refs. [25|,[26[ where a scaling of order n -1 logn 
was achieved by means of local adaptive strategies. In 
Ref. [25| it was also observed that if the channels 3> x 
are programmable 5 then a sub-shot-noise scaling of the 
accuracy is not allowed. Examples for which the no-go 
theorem of [2|| hold are provided by the class of clas- 
sical channels (i.e. quantum channels $ x which cannot 
propagate quantum information) and by the class of de- 
polarizing channels 30]: for these families of maps the 



4 For the Ramsey setup l^o"') corresponds to the maximally 
entangled state {\a)® N + \b)® N )/V2 of Fig. [TJ In the atomic 
literature, where \a) and \b) refer to different atomic levels, 

|i/>0 ) often is referred to as "spin-squeezed state" l42t l-Ufl or 
"Schrodinger-cat state" [441446H . It has been experimentally 
demonstrated using trapped ions [44H47H , Bose-Einstein conden- 
sates 48], macroscopic atomic ensembles [49|. and NMR |5C(|. 

5 The map Q x is called programmable if one can express it in 
terms of a constant interaction with an external ancillary system 
B whose initial state encodes the dependence on the parameter 
x, i.e. &x(-) = TrajVftT^ ® (')]V'} where V is a unitary trans- 
formation acting on B and S, rj x is the state of B, and Ttb{-} 
is the partial trace on B [5 j| . 
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best estimation accuracy of the parameter x will scale 
at most as n -1 / 2 . Recently it was also pointed out that 
for finite v, the two-step optimization approach adopted 
in the derivation of Eq. Q in general fails to provide 
the achievable bound [221 |52| . Instead one must adopt a 
min-max optimization scheme, minimizing the maximum 
RMSE, which always yields attainable bounds. It gives 
slightly worse performance than the two-step approach, 
but (apart from numerical prefactors) it maintains the 
same n^ 1 scaling (at least when the quantum channels 
<& x induce unitary transformations). In concluding, we 
also briefly mention that analogous scalings for spe- 
cific problems of quantum channel estimation have been 
obtained by using error measures which are substantially 
different from the RMSE adopted here (e.g. from those 
results one would not be able to infer the exact n~ l scal- 
ing for the RMSE which follows from the q-CR bound). 
In particular Refs. [33l - l37j consider the problem of es- 
timating unitary rotations in finite dimensional systems 
using a metric introduced by Holevo [TTj | . 



Applications in Quantum interferometry 

The prototypical example of a quantum interferomct- 
ric application of a quantum estimation procedure is pro- 
vided by the Mach-Zehnder interferometer. In this set- 
ting two input optical modes are merged at a first 50- 
50 beam splitter, propagate along two paths of different 
length accumulating an (unknown) relative phase shift tp, 
and are then merged at a second 50-50 beam splitter. The 
goal is to recover the value of tp by measuring the signals 
emerging from the interferometer, while employing a lim- 
ited amount of resources (i.e. by setting an upper limit N 
on either the maximum number or on the average num- 
ber of the photons entering the interferometer at each 
experimental run) . Indicating with a and b the annihila- 
tion operators associated with the two internal paths of 
the interferometer, the problem of recovering tp reduces 
to estimating a channel <5> x =ip which induces a unitary 
rotation e~ lxH with H = (a^a — tfb)/2 b eing the effec- 



tive system Hamiltonian, e.g. see Refs. 0,|!fl[53|. With 
this identification the two-step optimization strategy that 
brought us to Eq. ([5]) can be used to set a lower bound on 
the RMSE. First consider the situation in which generic 
POVM measurements are performed on v independent 
preparations of the interferometer. In this case Eq. ^ 
yields the following bound 



8tp v min 



|*> 2 y/V AH ' 



(6) 



where the minimization is performed over the set of in- 
put states |^) which satisfy the selected photon number 
constraint (i.e. the maximum number, or the average 
number constraint), and where AH is the associated en- 
ergy spread. Under both constraints, the optimal input 
1^) is provided by a state which at the level of the in- 
ternal modes of the interferometer, can be expressed as 



a "NOON" state USUI, a superposition of the form 
(| AT, 0) + |0, N))/y/2, in which N photons are propagat- 
ing al ong the first or the second optical path, e.g. see 
Refs. [5g]. NOON states are the formal analogue of the 
spin-squeezed states that achieve the Heisenberg bound 
in a Ramsey configuration (see footnote S] and Fig. [T]): 
in fact, also the NOON states exhibit a special sensitiv- 
ity with respect to the transformation which encodes the 
random variable tp (indeed NOON get transformed into 
output states (e-^ N / 2 \N, 0) + e lvN / 2 \0, N))/y/2 in which 
the phase tp result effectively multiplied by a factor N). 
With this choice Eq. © yields a lower bound of the form 



5tp v ^ 



/vN ' 



(7) 



which, for any given N , is achievable in the limit of 
large v, e.g. via maximum-likelihood estimation based 
on the photo-counting statistics at the outport ports of 
the interferometer. The bound (JJJ shows a A" 1 / 2 en- 
hancement with respect to more standard estimation ap- 
proaches where, for instance, the input ports of the in- 
terferometer are fed with coherent states of average pho- 
ton number A (these procedures show a SQL scaling 
5tp v = 1/y/v A in which, basically, all the vN photons 
contribute independently in the estimation process) . For 
this reason Eq. ([7J can be seen as the quantum optical 
counterpart of the Heisenberg bound of Eq. ((SJ). 

The attainability of the bound (J7J requires some extra 
considerations. First of all it assumes the ability of creat- 
ing NOON states. This is possible via some rather com- 
plicated optical schemes [57h61| which have only been 
implemented in highly refined, but post-selected exper- 
iments folTl66l ]. However states which possess a high fi- 
delity with the NOON states, also for large values of A, 
can be simply obtained by mixing a squeezed vacuum and 
a coherent state at a beam splitter |67H69j . Introducing 
such states at the input of a Mach-Zehnder, a scaling 
TV" 1 in the average photon number of photons employed 
in a given experimental run can be achieved To do 
so however one needs to employ the proper estimation 
process [70l |: a point which is not sufficiently stressed 
in literature where often suboptimal performances can 
be associated with a poor processing of the measure- 
ment outcomes. Analogous A^ 1 performances can also 
be achieved by employing different sources and/or by 
using estimator functions which are simpler than the 
maximum-likelihood approach. Typically these schemes 
are based on adaptive strategies where the parameter 
tp. is pushed toward an optimal working point which 
guarantees higher performances, e.g. two recent propos- 
als are Refs. [7l|, [72j while a list of older ones are in 
Refs. 4, 5]. Alternative schemes instead are based on se- 
quential strategies of Fig. [3] where a single photon pulse 
probes recursively the phase shift tp by passing multiple 
times through the delay line |24l - |26l . l38l |39| . 

The achievability of the bound (IT]) for finite values of v 
is also non trivial % E), HI M, HI, [M [zMI . Several nu- 
merical analysis supported the evidence that Heisenberg- 
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like TV -1 scaling should be achievable in the limit of large 
TV also for v = 1, where Eq. ([7]) would yield a scaling anal- 
ogous to the strong Heisenberg bound (see the previous 
section). In particular Refs. [5J, [7(J [3 [ill studied the 
asymptotic behavior of the confidence of the error prob- 
ability by adopting a Bayesian estimation strategy fill ]. 
A recent work by Hayashi [52| however appears to settle 
the problem by showing that while for v — 1 the bound 
of Eq. ((TJ is not exactly achievable, one could still reach 
an asymptotic N~ x scaling for the RMSE by adopting a 
min-max optimization approach [22j . 



Filtering protocols 

Instead of going through the trouble of creating the 
complex and fragile quantum states necessary for the 
quantum enhancements of quantum metrology, some pro- 
tocols have been proposed that use easy to create and 
robust classical states, and then filter (post-select) the 
high-resolution states at the measurement stage. The 
basic idea is to employ retrodiction [77], [zH : once a high- 
resolution quantum state has been detected at the out- 
put, one can interpret the whole experiment as having 
employed such a state since the input. This is a conse- 
quence of the fact that the wave-function collapse can be 
placed at an arbitrary time between probe preparation 
and measurement [z3, [79| . 

The filtering intrinsic in such protocols implies that 
part of the resources available at the onset are wasted: 
the system is sampled with many more resources than 
those actually employed for the parameter estimation. 
Moreover, since the system is sampled with classical 
states that present no quantum correlations, it is clear 
from the preceding sections that the SQL cannot be 
beaten: there is no increase in resolution over the opti- 
mal classical strategy that could employ all the resources 
(without filtering). For this reason, filtering protocols 
cannot be considered proper quantum metrology proto- 
cols, according to the definition we have given above. 
Nonetheless they can be extremely useful in the common 
case when the optimal classical strategy that employs all 
the resources that have sampled the system is impracti- 
cal. Moreover, in practical situations efficiency consider- 
ations rarely play a role, whereas robustness to noise is 
paramount. Classical states are, by definition, the most 
robust ones. In addition, there are situations where fil- 
tering methods achieve tasks that would be impossible 
with purely classical strategies, and post-selecting on the 
high-resolution quantum states is often very simple. 

In Ref. [78[ the theory of filtering protocols for phase 
estimation is developed, and necessary conditions are 
given to distinguish super-sensitivity (i.e. the error in 
the estimation is lower than allowed using classical re- 
sources) from super-resolution (i.e. the error in the esti- 
mation is lower than what would be allowed by a clas- 
sical procedure that uses only the resources that the fil- 
tering retains). Filtering protocols can achieve super- 



resolution (with an appropriate measuring strategy) but 
cannot achieve the super-sensitivity of quantum metrol- 
ogy proper. 

Beyond the Heisenberg bound: nonlinear estimation 
strategies 

Several Authors (80l496j have recently considered the 
possibility of using nonlinear effects to go beyond the 
TV -1 Heisenberg-like scalings in phase estimation prob- 
lems. These new regimes have been called "super- 
Heisenberg" scalings in Ref. [87[ , but a proper accounting 
of the resources shows that they are fully compatible with 
the analysis presented in the previous sections, see e.g. 
Ref. (97}. Ultimately the idea of these proposals is to 
consider settings where the unitary transformation that 
"writes" the unknown parameter x into the probing sig- 
nals, is characterized by many-body Hamiltonian gener- 
ators which are no longer extensive functions of the num- 
ber of probes employed in the estimation [83T - I92T ] or, for 
the optical implementations which yielded the inequal- 
ity ®j_ in the photon number operator of the input sig- 
nals |80l483l 193149a ]. Consequently, in these setups, the 
mapping ( e - l2: #)®™ which acts on the input states 

gets replaced by a transformations of the form e~ lxH ' 
which couples the probes non trivially. Accordingly the 
minimization of Eq. ^ is no longer forced to obey to 
the inequality ((7]). For instance, RMSE with scalings of 
the order oc N~ k can be obtained when using Hamil- 
tonians that involves fc-system interactions between the 
probes |85j . while oc 2"^ scaling can be achieved by 
introducing an exponentially large number of coupling 
terms [83J. Proposed implementations include scatter- 
ing in Bose condensates [13, [88| , Duffing nonlinearity in 
nano- mechanical resonators [931 ] . two-pass effective non- 
linearity with an atomic ensemble [§9| , Kerr-like nonlin- 
earities [80l - l82l l95l |96| , and nonlinear quantum atom- 
light interfaces [921 ] - 

Quantum metrology with noise 

The study of noisy quantum metrology is a special case 
of parameter estimation for channels, where the map <& x 
which describes the encoding of the unknown parame- 
ter x, contains also the description of the noise tam- 
pering with the process. Hence, even though very few 
general results are known, many results detailed in the 
previous sections can be used to characterize noise ef- 
fects. Non trivial examples are known where a Heisen- 
berg scaling can be retained even in the presence of 
noise 25]. On the other hand a simple application of 
the programmable criterion of Ref. [25[ shows that even 
a small amount of depolarizing noise is sufficient to ruin 
any sub-shot-noise performances one could reach by en- 
tangling the probing systems. Typical quantum metrol- 
ogy protocols are indeed extremely sensitive to noise. 
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For instance, the incoherent loss of a single photon in 
a NOON state transforms it into a statistical mixture 
(|JV- 1, 0)(N- 1, 0| + 10, N- 1)(0, N- 1|)/2 which is use- 
less for phase-sensing. Such extreme sensitivity to losses 
implies that this state cannot be used in any practical 
situation jfi^, [sH [§9[ . Is it then really possible to o utper - 
form classical strategies in practical phase sensing jlOOj ]? 

Surprisingly, it has been shown that asymp totically i t 
is possible to do so only by a constant factor [lOll Il02j : 
for any nonzero loss, for sufficiently high number of 
photons iV" the scaling of the optimal phase sensing is 
proportional to the scaling of the shot noise cx N^ 1 / 2 . 
While this means that quantum appr oaches are useful in 
highly controlled environments 102] (such as for grav- 
itational wave detection they only allow for very 
small enhancements in free-space target acq uisit ion [l02j j . 
Nonetheless, the shot noise can be beaten [l03| and the 
optimal states to do so in the presence of loss have been 
calculated numerically using various optimization tech- 
niques for fixed number of in put photons (l04L Il05j and 
for photon- number detection [l06j | . A post-selected proof 
of principle experiment that employs some of these op- 
timal states was recently performed (66|. Note that, for 
very lo w values of loss, NOON states retain their opti- 
mality }104ll06l |. and can be approximated by states that 
are easy to generate [67], H^- Also, a very simple proposal 
based on parametric downconversion which can be real- 
ized without post-selection was proposed in (72l |: it can 
achieve the Heisenberg bound for low loss and degrades 
gracefully with noise. In the case in which there is a large 
amount of loss after the sample has interacted with the 
light probe, even achieving the shot noise limit might be- 
come cumbersome (most of the photons that interacted 
with the phase shifter and contain phase information are 
lost). A simple strategy that amplifies the signal be- 
fore detection and can asymptoti cally achieve the shot 
noise was experimentally tested in [l07| . In contrast, the 
optimal states for the sequential (or multi-round) inter- 
ferometry [H, |39 | hav e not been obtained so far in the 
presence of noise [l00j |. However, as an alternative to the 
fragile two-mode states, some more robust single mode 
states were also analyzed, e. g. pu re Gaussian states in the 
presence of phase diffus ion Il08l| . mixed Gaussian states 
in the presence of lo ss fiol . or single mode variants of 
the two mode states [llfjj |. In contrast to the two mode 
case (where the phase is the relative one between the two 
modes) here the phase is measured relative to a strong 
classical signal (using heterodyne or homodyne measure- 
ments) or similar strategies. 

Historically, the first analysis of a quantum metrol- 
og y pro tocol [2j in the presence of noise was performed 
in where frequency measurements are analyzed in 

the presence of dephasing. Interestingly, a phase transi- 
tion was shown: for any nonzero value of dephasing, the 
maximally entangled state suddenly ceases to present any 
advantage over a classically correlate d sta te. (An anal- 
ogous result for magnetometry exists |112| .) We briefly 
present this result, as it is instructive of the subtleties 



encountered by quantum metrology when noise sources 
are considered. 

Consider the Ramsey setup of Fig.Q]when a dephasing 
process acts on the atomic probes. For the conventional 
setup the probability of finding the probe in the initial 
state after the application of the phase shift becomes 
p = (1 + e rt cosip)/2, where 7^0 measures the dephas- 
ing rate and t is the time elapsed from the state prepa- 
ration up to the measurement readout. Analogously, for 
the quantum setting we get p en t = (1 + e~ lNt cos Nip) /2, 
where the factor N in the exponent derives from the 
exponentially greater sensitivity of the entangled state 
to the dephasing. The corresponding RMSE (using the 

equations derived in Fig. [I]) is 6ip n — ( e ^ s 7n° s y ) 1 / 2 , and 



Sipn = ( 



,2~/Nt_ 



vN 2 sin 2 Nip 



-cos 2 Ntp \\/2 



) ' , respectively. For 7 > they 



are both minimized for t — (namely, we have to use 
a measurement procedure which is as fast as possible), 
and one recovers the typical N 1 / 2 enhancement of quan- 
tum metrology also in the presence of deph asin g. How- 
ever, if one wants to measure a frequency |2l. llll| ui, then 
the phase factor is ip = u)t. In this case, the errors of 
the separable and entangled procedure are respectively 
Su „ = { ^:- co if )Va and Su n = (i2^4p±y/\ 

n v nt 2 sin 2 uit I 't v sln 2 Nut ' 

For 7 = 0, we have S n uj — 1/ty/n and 8uj n = l/(tNy/U). 
The optimization over t says that our measurement must 
be as long as possible. Choosing the same t in the two 
cases, the entangled protocol is y/N more precise. In- 
terestingly, for 7 > the two situations cannot be any- 
more optimized on the same t. For example, if the fre- 
quency is ui — 7T7, it is easy to see that one has to choose 
t = (27) -1 to optimize the separable procedure and 
t = (2Nj)^ 1 to optimize the entangled one. With these 
choices, one obtains for the two cases 5uj n = y / e(2j) 2 /n 
and Suj n — W e(2j) 2 / v. We must assign equal resources 
to the two protocols, so we can repeat times the entangled 
procedure N times as it is N times faster. This recovers 
a factor 1/y/N also for the entangled procedu re, whose 
error then matches the separable one, Su n . In [ml fill 
it is shown that using non-maximally entangled states 
one can beat the separable procedure by ~ 40%, which 
was proven to be the ultimate gain. 

Summarizing, the estimation of the phase ip is unaf- 
fected by the presence of dephasing noise (as long as 
the measurement time t can be chosen short enough): 
the full N^ 1 / 2 enhancement of quantum metrology is re- 
tained. On the contrary, the estimation of the frequency 
u = ip/t is strongly affected by the same noise: basically 
all quantum enhancement is lost for any value of 7. This 
is connected to the fragility of the entangled resource in 
the presence of noise: the entangled procedure must be 
performed N times faster, which reduces the precision 
in the estimation of uj by a factor N, thus achieving the 
same final yield as the separable protocol. Note, how- 
ever, that if the observation time t in atomic clocks is 
constrained by experimental issues (typically, by the lo- 
cal oscillator fluctuations}, then entanglement all ows a 
sub-shot noise scaling 114 which scales as N~ 2 / 3 [l!5| |. 
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A generalizat ion of the frequency measurement theory 
was given in |l!6j |. which covers many noisy estimation 
measurements when the measurement duration time and 
production rate of the probes are considered as relevant 
resources. 
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